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Abstract 

We propose an alternative method to study the saddle point equation in the strong coupling limit for the Wilson 
loop in Af = 2 D = 4 super Yang-Mills with an SU(N) gauge group and 2N hypermultiplets. This method is 
based on an approximation of the integral equation kernel which allows to solve the simplified problem exactly. 
To determine the accuracy of this approximation, we compare our results to those obtained recently by Passerini 
and Zarembo. Although less precise, this simpler approach provides an explicit expression for the density of 
eigenvalues that is used to derive the planar free energy. 
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1 Introduction 



One way to obtain non-perturbative results in supersymmetric gauge theories is to use the localization technique to 
rewrite some of the observables, such as the partition function or the Wilson loops, as multivariable integrals [1]. 
These integrals are reminiscent of the eigenvalue integrals that appear in matrix models [2]. The number iV of 
integration variables is associated to the size of a gauge group and the large N expansion interpreted as a 't Hooft 
expansion. To compute the integrals order by order in this l/N 2 topological expansion, various techniques can be 
borrowed from matrix models. In particular, the planar order can be obtained using a saddle point approximation: 
the integration variables, or 'eigenvalues', condense on a continuous support and the saddle point equation takes 
the form of an integral equation for their density. This approach has proved to be very efficient to study three 
dimensional supersymmetric gauge theories, such as Chern-Simons or ABJM (see [3] and references therein). 

In four dimensions, the Wilson loop and free energy of the N = 4 super Yang-Mills theory are described by the 
standard Hermitian matrix model with Gaussian potential. This was first conjectured in [4, 5] after a perturbative 
analysis, and later proved by Pestun using localization [1]. In this simple case, the eigenvalue integrals can be 
performed exactly, leading to an expression of the Wilson loop valid at all order in N and for all values of the 't 
Hooft coupling A = Ng 2 , where g is the Yang-Mills coupling. This gauge theory is known to be dual to the type 
IIB string theory in AdS$ x S5 background [6]. At large A, the Wilson loop is characterized by an exponential 
behavior in \/X at the planar order, which agrees with the supergravity result [4]. The subleading orders were also 
successfully compared to the dual string amplitude in [5]. 

The N = 2 Ad super Yang-Mills theory with an SU(N) gauge group and 2N fundamental hypermultiplets, also 
referred as superconformal QCD (SCQCD), is more intriguing. Although several proposal were made [7-10], the 
dual string background has yet to be identified. On the gauge theory side, localization can still be used to express 
the partition function and Wilson loop as integrals over real eigenvalues. But the Van der Monde determinant of 
the Hermitian matrix model is now replaced by a product involving G-Barnes functions arising from the one-loop 
contribution [l]. 1 This model was first proposed in [12], where a scaling argument has been employed to predict 
a very unusual power growth of the planar Wilson loop at large A. This fact was later confirmed by Passerini and 
Zarembo who used the saddle point technique to derive the expression of the Wilson loop at strong coupling, 

W{\) = a( X \ , A~9.47xl0~ 5 . (1.1) 

VviogAy 

We briefly summarize their computation below. 

In the saddle point approximation, the planar Wilson loop is written as an integral over the density of eigenvalues, 

W(X) = J e 2nx p(x)dx, (1.2) 

which satisfies a singular integral equation, the saddle point equation, 

I p(y)dy _ f k^ x _y^p^dy = ^L x -k(x). (1.3) 

The RHS of this equation is the derivative of the matrix potential, but it will simply be referred here as the potential 
of the integral equation. The kernel of the equation (1.3) is a linear combination of the usual Cauchy kernel and 
the function k. In the M = 4 theory, the function k doesn't appear and we recover the saddle point equation of the 
Gaussian Hermitian matrix model. The solution is then given by the Wigner semi-circle law, and the endpoints of 
the compact support V scale as y/X. This square root dependence leads to the Wilson loop exponential behavior at 



'The instanton corrections are expected to give exponentially decreasing corrections at large N [1 1] and will be neglected in this note. 
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large A. In SCQCD, the function k is the logarithmic derivative of the G-Barnes functions that replace the Van der 
Monde determinant, it can be expressed as a sum of two digamma functions, 

k(x) = x(ip(ix) + ip(-ix) + 27), = -7. (1.4) 

In the weak coupling limit A — > 0, the function k can be neglected and we recover the N = 4 solution. Conversely, 
at strong coupling this function dominates the Cauchy kernel, leading to an unusual behavior of the density. To in- 
vestigate the saddle point equation (1.3), the authors of [11] employed the Wiener-Hopf method, a refinement of the 
Fourier transform for integral equations with compact support. They derived an approximate density in the region 
x > which was expressed in the Fourier space. From the unit norm condition, and using the parity of the density, 
they deduced a logarithmic dependence of the support endpoints in A. This exotic behavior of the support gives rise 
to the algebraic expression (1.1) of the Wilson loop in A and log A at strong coupling. Their method and results were 
consistently checked by numerical simulations. 



A possible approach to uncover the string background of SCQCD is to consider M = 2 super Yang-Mills 
with SU(N) x SU(N) gauge group and bifundamental hypermultiplets. This theory depends on two 't Hooft 
couplings, one per gauge group, and the ratio of the two deforms continuously a Z2 orbifold of N = 4 into the 
N = 2 SCQCD. The dual string background of this theory is known to be AdS^ x S5/Z2, and some insight into the 
SCQCD dual background was obtained in [8] from this (singular) limit of the two couplings ratio. After localization, 
the SU(N) x SU(N) gauge theory is described by a quiver matrix model: the two Wilson loops, associated to each 
gauge group, are coupled trough their saddle point equation [12, 13]. Although remarkable, the calculation presented 
in [1 1] may be difficult to extend to the quiver matrix model. We propose here a simpler method to study the saddle 
point equation (1.3) at strong coupling that still retains the qualitative aspects of the problem. Keeping in mind a 
later application to the quiver problem, we will systematically compare our results with those derived in [1 1] in order 
to determine the accuracy of the approximation. 

It can be shown by a rescaling of the eigenvalues that the large A limit amount to send the spectral parameters 
x, y to infinity [12]. We suggest to truncate the function fc's large x expansion in the kernel and potential of (1.3), 

00 „ 

k(x) = 2xlog|x| +2 7 x + ^(-l) fc -^x~( 2fc+1 \ x^±oo, (1.5) 

where denote the Bernouilli numbers. The replacement of the kernel is justified for a large distance interaction 
between eigenvalues, i.e. when \x — y\ is large. When the eigenvalues are close to each others, the kernel k is 
negligible compared to the Cauchy kernel, and the effect of the replacement should remain small. The approximation 
of the potential for \x\ large pertains to the Wilson loop computation: the dominant term is given by the largest 
eigenvalue which tends to infinity at strong coupling. Although convenient, this second approximation is not strictly 
necessary, and the solution for the full potential will be given in the concluding section. 

We restrict ourselves to the first orders of the expansion (1.5), for which the integral equation can be cast into 
the form 2 

/ \-a / k sc (x - y)p(y)dy = fix + ak sc (x), k sc (x) = x log \x\ - x. (1.6) 

Jr x-y Jr 

This simplified problem can be solved exactly, providing an explicit expression for the density. Although less precise, 
this expression may appear more convenient than the one found in [11]. It will be used to compute the Wilson loop 
with the formula (1.2). It can also be employed to derive the planar free energy which is the effective action taken at 
the saddle point, 

= -S eS [p] = - J p(x)dx (~Y x2 - 2L(x)^j + p(x)dx p(y)dy (log \x - y\ - L(x - y)) (1.7) 



2 We use here the invariance of the equation under shifts of the kernel k by a linear function which is a consequence of the density parity 
(see section 2 below). 
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where L is the logarithm of the G-Barnes functions, L' = k and L(0) = 0. 

To have a hint on the effect of higher terms in the kernel truncation, we consider two different orders of approx- 
imation. At the first order, we take 

g 2 

k{x)~2k tR {x) + 2(i + \)x, a 1 = -2, fa = — . (1.8) 

As arefined approximation, we also take into account the 1/x term of the expansion (1.5) which is of the same order 
than the Cauchy kernel, 

k(x)~2k sc {x) + 2(j + l)x + ^, a 2 = -^, /3 2 = -^. (1.9) 

X OA 

The effect of the subleading term Bijx in the potential, i.e. in the RHS of (1.3), will be neglected here as it remains 
subdominant compared to all other members. In the following, we use the notation a = —a 2 (with a > 0) and 
(3 = 4ir 2 a 2 /\. The two approximations (1.8) and (1.9) will be referred to by the indices 1 and 2 respectively. 

In this setting, we keep the linear term (3x that comes from the Gaussian matrix potential and confines the 
eigenvalues. Consequently, the support of the density remains compact and the Wilson loop is finite. This method 
keeps all the general characteristics of the solution observed in [11]. In particular, we recover the logarithmic 
dependence in A of the support endpoints. The Wilson loop we computed exhibits a similar behavior than (1.1), with 
the approximate value 3.056 of the cubic exponent (using the refined approximation (1.9)). It is argued that the exact 
exponent appears as a limit where all the terms of the expansion (1.5) are considered, and that this can be encoded 
in an effective value for the variable a, = tt/2. 

From the solution of the approximate integral equation (1.6), we deduced the first orders of the planar free energy 
at large A, 

16vr 2 16vr 4 (logA) 2 /logA^ 



^ (A = oo)--^ + ^ r ^ + O(- ?r j (1.10, 

where Jo (A = oo) is an integration constant we determined. This is one of the main results of this note, together 
with the Wilson loop expression (2.32). The form of the A-dependence in (1.10) does not involve a, it is believed to 
be identical to the exact result that would be obtained from the full density solving the initial problem (1.3). 

2 Solution of the integral equation 

To solve the simplified integral equation (1.6), we notice that, like the original equation (1.3), it is invariant under 
the sign inversion of the eigenvalues x — > — x. As a result, we look for an even density with a symmetric support 
T. Since the eigenvalues are real, this support lies on the real axis, and we concentrate on a single cut solution for 
which T = [— M, M]. The density p has to vanish at the endpoints of its support, and is normalized to one. Due to 
the density parity, the first moment is vanishing which implements the SU (N) constraint. 

The first step to solve (1.6) is to transform the Cauchy kernel into the k sc kernel. It is done using an integration 
by parts, 

I p(y)dy = I _ y)p „ {y)dy _ pl{M) [hc{x + M) + hc{x _ M)] (21) 
Jt x ~ U Jr 

We will see below that p'(M) is actually infinite, and this identity only makes sense after a proper regularization. 
We will elaborate on this subtle point later, and for now we formally use (2.1) to rewrite the saddle point equation as 
an integral equation with kernel k sc over the auxiliary density ip = p" + ap, 

J k sc (x - y)<p(y)dy = fix + ak sc (x) + p'{M) [k sc (x + M) + k sc (x - M)] . (2.2) 
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In doing so, the potential acquires a boundary term. Since ip is even, this equation is equivalent to its derivative, 

J log \x - y\(p(y)dy = (3 + a log \x\ + p'(M) log (M 2 - x 2 ) . (2.3) 

This weakly singular integral equation is linear in 92, and the full solution is a sum of three components associ- 
ated to the three members of the potential. For each of these solutions, that we denote tpi, we have to solve an 
inhomogeneous differential equation to recover the associated density pi, 

p'l (x) + a Pi {x) = tpi(x). (2.4) 

The integration constants are fixed by imposing the parity of the density and the vanishing condition at the support 
endpoints. Below, we discuss each solution pi individually and it is convenient to fix the integration constants 
independently on each p, L . 



2.1 Solution for the boundary term 

We first look at the solution for the boundary term, i.e. the last term in the RHS of (2.3). This solution pj is simply 
a sum of two 5 functions supported on the endpoints x = ±M. The differential equation (2.4) for the density can be 
solved using a Fourier transform, and we get 

p/ ( x ) = -P^Ml L-a\x+M\ + e -a\x-M\ _ 2e ~aM CQshax \ . (2 5) 

This density is vanishing on the whole support V and cannot contribute to any regular integral. Consequently, it has 
a zero norm and gives no Wilson loop contribution. As we will see in the subsection 2.3, the main role of pj is to 
cancel the divergent terms that appear due to the integration by parts we employ. The derivative of pj is supported 
on the two endpoints of T where it takes the values 

f/j(±M) = ±\p'{M). (2.6) 

The unknown quantity p'(M) in the expression (2.5) of pj will thus be determined once we know the two other 
components pjj and pjjj of the density p, 

l -p'(M) = p' II {M)+p' III {M). (2.7) 



2.2 Solution at (3 = 

The density ipn solves the integral equation (2.3) for the logarithmic term in the potential. It is given by a delta 
function, now supported at x = 0, which leads to the density 

, . a sinha(M — Ixl) , n ns 

"'<*> = 2 coLm <2 - 8) 

after solving the differential equation. At infinite coupling /3 = 0, pjj is the full solution, up to the boundary term 
pj. As already observed in [11], the unit norm condition on pjj leads to a non-compact support, 



/ p n (x)dx = 1 — — — — ■ = 1 M ->• oo. 

J T ' cosh(aM) 



(2.9) 



3 An integral equation with the kernel k sc and a constant potential cannot have an even solution, so that the integration constant is automat- 
ically zero. 
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This is due to the vanishing of the Gaussian potential that confines the eigenvalues. In the next subsection, we 
introduce the density pm associated to this confining potential, i.e. to the constant term in the RHS of Equ. (2.3). 
This density pm contributes to the norm of p, and M remains finite. 

Since at j3 = the support F = R is non-compact, the initial equation (1.3) with complete kernel can be solved 
by Fourier transform [11]. It gives us the opportunity to compare our solution for the approximated kernel 

p n ( x ) = | e - a N at M -> oo, (2.10) 

to the exact solution of the problem, 

Poo(x) = - ) = e-W/ 2 + O (e-^) , (2.11) 

2cosh(7rx/2) V / 

in the strong coupling limit \x\ — > oo. We note that the coefficients inside the leading exponentials are different, 

tt/2 -1,571, ai = y/2~ 1,414, a 2 = 2y^3/5 ~ 1, 549, (2.12) 

but the numerical discrepancy is rather small (1.4% for the finest approximation (1.9)). This discrepancy can be 
easily explained in the Fourier transform framework. Our strong coupling approximations of the kernel (1.8) and 
(1.9) amount to expand the Fourier transform k(uj) of k at small ui. Then, the two Fourier transformed densities agree 
at the first orders poc{w) — Pii(w) but this is not enough to ensure the correct asymptotic in the x-space. However, 
given the proximity of the numerical value in (2.12), we expect a good approximation for the cubic exponent in the 
Wilson loop expression (1.1), with an error of only a few percents. 

We further notice that the coefficient a in (2.12) seems to converge toward the effective value = tt/2 as we 
refine the approximation for the kernel. To justify this fact, we may consider the next orders in the expansion (1.5) 
of the function k. After integration by parts, it leads to define an auxiliary density involving higher derivatives in 
p. Neglecting boundary terms, this density tp still satisfies the equation (2.3) with /c sc kernel and a similar potential. 
When taking into account the whole expansion (1.5) for k, the auxiliary density writes 

<p = P " - 2 P - ± ^7 2 rr/ 2m+2) - (2 - i3) 

m=0 v ' v ; 

The coefficient a is related to the homogeneous solution of the differential equation and can be obtained as a root of 
the characteristic equation. The characteristic equation for the infinite differential system (2.13) takes the form of an 
infinite series. Remarkably, the summation can be performed and we are led to the very simple equation, 

2sin 2 (a/2) = 1. (2.14) 

The solutions are aoo = tt/2 modulo it. Since the density should vanish at \x\ — > oo, we have to select only the 
decreasing exponentials, and we recover the dominant contribution of (2.11), as well as all the exponentially small 
corrections. 

2.3 Complete solution 

We now analyze the last part of the density, pm, which is associated to the constant term /3 in the potential of (2.3). 
The integral equation for the auxiliary density (pjjj does not depend on a, and we first take a = 0. Then, the integral 
equation (1.6) on p reduces to the N = 4 saddle point equation which is solved by the Wigner semi-circle law, 

p{x) = -R(x), R(x) = VM 2 - x 2 . (2.15) 

7T 
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At a = 0, ipiu is equal to the second derivative of the Wigner density, up to a boundary term <J~^\, 

^ /(X) = -J^ + ^ (X) - (2 " 16) 

Pluging this ansatz into the integral equation (2.3), we indeed find a solution for ip H ' I as a sum of two delta functions 
supported on the endpoints, 

= Jw) {6{x + M) + 6{x - M)) - (2 - 17) 

This expression, obtained after integration by parts, contains the infinite term 1/R(M). There are several ways to 
regularize this computation, the simplest is to replace the Wigner solution R(x) by the function 

R e (x) = V(M + e) 2 -x 2 . (2.18) 

The corresponding density provides the correct solution up to O(e) terms, and R e is non-vanishing at x = ±M. To 
avoid unnecessary complicated expressions, we will not explicitly write down this regularization here, and we keep 
the notation 1 / R{M) to imply 0(1/ y/e) divergent terms that have to be canceled. This cancellation will be achieved 
using the density pi, which shows that the infinite terms appearing here are just an artifact of the integration by parts. 

Returning to a ^ 0, we have to solve the differential equation (2.4) on pm with the second member (pm given 
by Equ. (2.16). We get 



/3M 2 
Piii{x) = -- 



[ x smha(x-y) _ cosh ax f M sinho(M-y) 1 {B ) () (2 , g) 



r-M 

j sum uyx ■ — y) ^ uusnux j 

ira 

and a boundary term which is similar to the expression (2.5) of pj, 



pi ^ x) = -^iw) ( e ~ a]x+Ml + e ~ a]x ~ Ml ~ 2e ~ aM coshax ) ■ (2 - 20) 

We used here the homogeneous solution to impose the vanishing of the density at x = ±M. The expression (2.19) 
for pm is even and reduces to the Wigner solution at a = 0. 

Once pi i and pm are known, we can determine the constant p'(M) involved in the expression of pi using the 
relation (2.7). The derivative of pm at the endpoints contains a 1/R(M) singularity, 4 

where l\ denotes the Bessel function of the first kind. From this result and the derivative of pn at x = M, 

( B) 

we deduce the expression of pi. As required, the divergent terms in pi and p HI cancel each others, and we are 
led with a finite expression for the density. This final density can be written as a sum of three components, p = 



4. 



This calculation is a bit delicate and we have to use the formal formula 

M 



f cosh ay , cosh(aM) m , , 



6 



Pi + Pi I + Pi 1 1 with 



Pi(x) 
Pn(x) 
Piii(x) 



a-PMhjaM) f- a \ x+M \ a \ x _ M \ _ ^ aM \ 
2 cosh aM V J 
a sinh a{M — \x\) 



2 cosh aM 
(3M 2 \ f x sinha(x - y) 



(2.24) 



ira 



f x sinh a 

Jo R( 



R{yf 



dy 



cosh aM 



f-M 

Jo 



cosh ax f sinh a (M — y) 



Riyf 



dy 



We still have to determine the dependence in A of the support T, which is usually achieved by imposing the unit 
norm condition on p. The norm of p have contributions only from pjj (Equ. (2.9)) and pm, 

j r Pin(x)d X = acQshaM . (2-25) 

Combining these two results, we deduce the relation 

aMh(aM) = -4r. (2.26) 

When A is infinite, M should also be infinite, as seen in (2.9). Using the asymptotic of the Bessel function, the 
previous relation simplifies into 

CVMe aM = A M (A) ~ - log A - — log log A + • • • (2.27) 

a 2a 

where C is a constant factor. We recover the logarithmic dependence of the endpoint M(A) found in [11] and 
specific to SCQCD. More precisely, the relation between M and A is identical to the Equ. (4.7) of [11], provided we 
use the effective value a^ = 7r/2. The constant C should be compared to the value 14.60 obtained in [11], and to 
the outcome of the numerical simulation Cint — 15.02 (appendix B of [11]) 

C = (2tt) 3 / 2 ^, d ~ 18.73, C 2 ~ 19.60, ~ 19.74, (2.28) 

where the index oo refers to the use of the effective value a^. We note that our method gives the correct order of 
magnitude for C, but does not provide a very accurate approximation for this constant. 

2.4 Wilson loop 

The Wilson loop is computed using the formula (1.2) with the solution (2.24) of the density. There are two contri- 
butions, 

a 2 /cosh2^M 

WllW = A^a~ 2 \ cosh aM ~\ 

V x (2.29) 

,„ PM ( T tn cosh27rM w 

These expressions are exact in the approximation (1.6) of the saddle point equation. Interestingly, the dominant 
term at strong coupling, Win, exhibits a compact expression as a Bessel function of M. This is reminiscent of 
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M = 4 super Yang-Mills [4,5], and one could speculate about a Laguerre polynomial formulation for the l/N 
corrections [5]. At the first order, the Bessel functions can be approximated by their asymptotic, 5 



W(\)^r^-e 2 * M , with r= (2 " )3/V 2 ( ^-^ ) . (2.31) 
A 4tt z — or 

The dependence of the Wilson loop in M remains exponential, and the power growth of W{\) is due to the unusual 
logarithmic behavior (2.27) of endpoint M(A), 

w W *a(*)*-\ A= w'*-»'-fyg-^, (2 . 32) 

V V log A / 47r z - a A 

Using the refined approximation (1.9), we get the exponent 3.056, which differs only by 2% from the cubic exponent 
found by Passerini and Zarembo. This cubic exponent arises at the effective value = tt/2 which takes into 
account the neglected higher order terms in the kernel approximation. Our approximation becomes less accurate for 
the constant factor A, 

A x ~ 4.46 x 1(T 6 , A 2 ~ 1.44 x 1(T 5 , A^ ~ 1.71 x 10" 5 , (2.33) 
which differs roughly by a factor five from the value A ~ 9.47 x 1CT 5 derived in [11]. 

2.5 Free energy 

To find the planar free energy, we should inject the expression of the density into the effective action (1.7). But it is 
actually easier to first consider the A-derivative which is simply the second moment of the density, 

9\MP] = J^x 2 p{x)dx. (2.34) 

This relation can be obtained by taking the derivative of Fq [p] seen as a functional of p, and then use the solution 
p = p* of the saddle point equation. It can also be checked explicitly that 



d_ 

dX 



(T [p*}) = (d x T )[p*] -2 j^d x p*(x)V eS (x)dx, (2.35) 



where 



V efl (x) = -^-x 2 - L(x) - p{y)dy [log \x-y\- L(x - y)] (2.36) 

is the effective potential associated to the integral equation (1.3). Note that the A-dependence of the support gives 
no contribution since p(M) = 0. When p satisfies the saddle point equation, the effective potential is a constant and 
can be removed from the integration in the RHS of (2.35). This term becomes proportional to the derivative of the 
density norm, which vanishes since p has norm one. 

The second moment of the density remains finite at strong coupling, and the leading order of the free energy 
A-derivative (2.34) can be computed using the infinite density profile p^ given in (2.11), 

dx^o = -j^r / x 2 poo(x)dx + o ( -^J = +o (-^ J ■ ( 2 - 37 ) 



5 Here the coefficient r differs roughly by a factor 2 from the result r ~ 2.18 derived in [1 1], 



n ~ 1.11, r 2 ~ 1.29, r x ~ 1.32 . (2.30) 
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This expression can be integrated to obtain Fq, 



jf = jr (A = oo) - ^ + o QY (2.38) 

where the integration constant is the free energy at infinite coupling. To find this constant, we notice that the insertion 
of the effective potential in the expression (1.7) of the planar free energy leads to a huge simplification, 

Jb(A) = j p *{ x )\og\x\dx- j x 2 p*(x)dx. (2.39) 

At infinite coupling, the second term vanishes and we can use again the density p^ to compute the integration 
constant, 

Jo(A = oo) = jf Poo (s)log|s|dx = 2 log (y^) ~ -°- 7832 - ( 2 - 4 °) 

The first and second orders of the planar free energy at strong coupling, resp. (2.40) and (2.38), were computed 
exactly and no approximation of the integral equation (1.3) has been used. To go beyond these first orders, we need 
to take into account the correction to the density due to the confining potential of the integral equation. This is 
exactly the role played by pm for the strong coupling approximation of the integral equation. Thus, the density 
(2.24) determined previously should give us a hint on the form of the subleading term for the free energy, but at the 
price of a distortion on the numerical coefficients. In particular, p^ will be replaced by pjj. 6 

The relations (2.34) and (2.39) are exact provided that p solves the initial saddle point equation (1.3). When 
employed with the density solving the equation with approximate kernel (1.6), it assumes a similar strong coupling 
approximation for the expression (1.7) of ^b[p]- The calculation of the second moment for the density (2.24) is 
relatively straightforward, and using the normalization condition (2.26) the free energy derivative simplifies into 



2 



16tt 

^A 2 



Hi- 2 1 32vr 4 (logA) 2 



1 



2ir 2 a 2 M' 



a 2 A 2 a 2 A 3 



+ ■■■. (2.41) 



A 

Using the relation (2.27) between M and A, the expression (2.41) can be integrated to obtain To, 

Jo ^ W = oo) - ^\ + ^(1 + aM)e~ 2aM 

° C2 42 N l 

. 16vr 2 l 167r 4 (logA) 2 „{\og\^ 

As expected, the first and second order terms have the same form than the exact result, with different numerical 
coefficients. The next order contains a logarithmic non-analyticity that would be interesting to interpret in the gauge 
theory context. To test the distortion of the numerical coefficients, we computed the integration constant for the 
approximate density pjj, 

J 7 o(X = oo)= / pn(x) log \x\dx = — 7 — loga2 c± — 1.015, (2.43) 

the numerical discrepancy between the exact result and the approximation is about 30%. The subleading coefficients 
differ by 17% in the approximation (1.9). 



6 Computing the free energy, the potential approximation at large \x\ is not fully justified. However, we shall see in the next section that it 
does not modify the form of the subleading term. 
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3 Discussion 



In this note, we discussed an alternative method to study the saddle point equation for the Wilson loop in SCQCD. 
We used a strong coupling approximation for the kernel and potential of the saddle point equation, and derived 
an explicit solution for the density of the simplified problem. From this solution, we computed the Wilson loop 
dominant contribution at strong coupling and compared it to the expression derived in [11]. We found a good 
numerical approximation for the exponents characterizing the dependence in A and log A. The exact exponents can 
be obtained from an effective parameter that takes into account the effect of neglected terms in the kernel strong 
coupling approximation. The simplified problem exhibit the same qualitative behavior than the original one, and our 
approach should be suitable to simplify the quiver matrix models when at least one of the gauge couplings is strong. 
A similar technique can also be employed to study the expressions of the 't Hooft loops recently derived for N = 2 
gauge theories from localization [14]. 

We computed the leading orders of the free energy, the subdominant term exhibits an intriguing 1/A behavior 
that differs from the logarithmic result found in N = 4 super Yang-Mills theory, 

4^ =4) = ilogA-i-log(47r). (3.1) 

Reintroducing the variables N and g, and taking into account the 1/N 2 factor in front of the planar free energy, 
the SCQCD free energy behaves at strong coupling as T ~ 1/iV 2 + l/g 2 N 3 (we omitted the constants). This 
information could help to select the correct background of the dual string theory. 

Until now, we mostly focused on the kernel approximation. To investigate the approximation of the potential we 
consider 

f + a f ksc(x _ y )p(y)dy = (3x + % (k(x) - 2( 7 + l)x) (3.2) 

JY x ~ V JT * 

where we kept the full potential of (1.3). Since the integral equation is linear in the density, we only need to introduce 
an additional contribution pjy to the total density p to take into account the higher order terms. Using the asymptotic 
expansion (1.5), it is not hard to see that 

PMx) = "Y EtTTT^T E « m (-l) fc - m <5( 2fc - 2m >(*) + \f(a) P n(x). (3.3) 

k=0 ^ '' m=0 

with 

00 T3 ( 1 \k 2 

f(a) = V B2k+2 a 2k + 2 - ° 2 (3 4) 

h k + l ^ "2sin 2 (a/2) (J - 4) 

solves the problem. This density contributes to the normalization of p, which modifies the constant C in the relation 
(2.27) giving M(A), 

C = 4(2^) 3/2 a" 3/2 sm 2 (a/2), d ~ 15.81, C 2 ~ 15.98, = 16. (3.5) 

These figures are closer to the value C = 14.60 found in [11]. The contribution Wiy to the Wilson loop remains 
subdominant, but the modification of the relation M(A) affects the prefactor A in the expression of W(X), 

A! ~ 9.48 x 1(T 6 , A 2 ~ 3.30 x 10~ 5 , Ax> ^ 3.95 x 10" 5 , (3.6) 

which now only differs roughly by a factor two from the numerical value given in [11]. The density pjy also 
contributes to the derivative of the free energy, and we get in the approximation (1.9), 



2 



AT 87T 



1 



4^ 2 M 2 



8tt 2 327r 4 (logA) 2 



A 2 a 2 A3 



+ • • • (3.7) 
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We notice that the dependence in A remains the same when we consider higher corrections for the potential, but the 
overall constant factors get modified. In this way, we reach a better agreement when comparing to the results of [1 1]. 
But these numerical factors remain rather sensitive to the approximation, and we cannot rely on our approach for 
their derivation. On the other hand, the A-dependence is quite strong and can be trusted. 

So far, we mainly concentrate on the technical aspects of the problem, and these results still have to be understood 
within the gauge theory framework. In particular, it would be interesting to see how the logarithmic non-analyticities 
appear at strong coupling. In this context, we cannot employ the kernel approximation, and should study the original 
equation (1.3), which solution interpolates between strong and weak coupling regimes. The kernel function k has 
poles on the imaginary axis at x G iL. To each pole is associated a copy of the eigenvalue density support where the 
resolvent of the corresponding Riemann-Hilbert problem has a branch cut. In the weak coupling limit, x is small and 
the branch cuts associated to the imaginary poles are send to infinity. Only the branch cut T lying on the real axis 
remains. On the other hand, at strong coupling these branch cuts condense on the real axis, leading to the logarithmic 
behavior for the kernel. One may wonder if something similar happen within the gauge theory. 
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